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ABSTRACT 

This  paper  evaluates  five  maneuvering  strategies  for  an 
evader  in  a  two  dimensional  continuous  evasive  game.   An 
evader's  movement,  in  two  dimensions  at  constant  speed,  is 
simulated  by  choosing  courses  that  are  independent  and  dis- 
tributed according  to  one  of  the  five  course  change  rules. 
The  times  between  course  changes  are  independent  exponential 
random  variables.   An  improvement  over  a  previously  estab- 
lished least  upper  bound  on  evader  survivability  is  pre- 
sented.  The  bulk  of  this  paper  is  the  presentation  of 
evader  survivability  as  a  result  of  the  maneuvering  strategy 
he  employs . 
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I.   INTRODUCTION 

A  classic  military  problem  with  current  application  is: 
how  to  maneuver  a  mobile  target  in  order  to  prevent  success- 
ful prediction  of  its  position.   For  example  how  should  an 
aircraft  carrier  maneuver  to  maximize  its  survivability 
against  an  ICBM  attack? 

Of  primary  importance  to  this  problem  is  the  time  lag 
between  the  decision  to  fire  at  a  target,  made  at  time  t, 
and  the  warhead  detonation  at  time  t+TL.   This  lag  is  called 
time- late  and  denoted  as  TL.   The  decision  to  initiate  such 
an  action  is  based  in  part  upon  knowing  the  target's  posi- 
tion at  time  t  while  the  effectiveness  of  that  action  is 
dependent  upon  the  error  in  predicting  that  position  at  t+TL. 

A.   ICBM  VERSUS  AIRCRAFT  CARRIER 

Consider  the  ICMB  versus  aircraft  carrier  problem.   As- 
sume that  the  carrier  operates  in  mid-ocean,  possesses  great 
endurance,  can  make  sharp  turns  and  has  as  its  only  kinematic 
restriction  a  constant  speed,  v.   The  carrier  is  referred  to 
as  the  evader,  E.   The  attacker,  P,  continuously  observes  E 
from  a  nearby  unarmed  trawler.   P's  weapon  is  a  land  based 
ICBM  which  has  perfect  accuracy,  produces  a  lethal  area  A 
(this  lethal  area  can  have  any  shape)  but  does  not  have  a 
mid-course  guidance  nor  homing  capability.   Although  P  knows 


E's  position  past  and  present,  he  also  knov/s  that  his  ICBM 
will  have  a  one  hour  time-late. 

P's  problem  is  to  predict  E's  position  at  time  t+TL  in 
order  to  initiate  an  attack  at  time  t.   E's  problem,  since 
he  will  not  know  of  an  attack  \mtil  he  observes  the  deto- 
nation, is  to  maneuver  in  such  a  way  to  confound  P's  predic- 
tion.    E's  movement  should  be  random  because  any  nonrandom 
movement  would  be  vulnerable  to  extrapolation. 

B.   AN  UPPER  BOUND  ON  EVADER  SURVIVABILITY 

A  lower  bound  on  P's  success,  p,  ,  and  consequently  an 
upper  bound  on  E's  survivability,  l~Pvf  is  known.   P  knows 
E's  position  at  t+TL  will  be  interior  to  or  on  a  circle  of 
radius  vTL  centered  on  E '  s  position  at  t  (which  is  known). 

The  enemy  can  then  gaurantee  a  kill  probability  of  at  least 
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A/tt(vTL)   by  choosing,  a  lethal  area  randomly  within  the 

circle  of  uncertainty  (a  wedge  of  random  orientation  will 

do)  [Ref .  4] .   Then  the  maximum  survivability  E  could  attain 

2 
would  be  1-A/7t(vTL)  . 

This  upper  bound  on  E's  survivability  is  a  function  only 

of  P's  strength,  that  is  the  magnitude  of  A  (the  value 

2 

tt(vTL)   is  constant).   Figure  1  is  a  graph  of  that  survival 

function  versus  P's  probability  of  success. 


Time-late  is  a  summand  of  many  factors  in  the  command  and 
control  problem  as  well  as  the  missle  flight  time.   For  the 
purpose  used  here  it  is  sufficient  to  treat  time-late  in 
total  and  not  consider  its  decomposition. 
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If  E  knows  the  time  at  which  an  attack  is  initiated  then 

the  problem  is  trivial. 


0.5         1.0   pj^ 
Survivability  Versus  Attacker  Strength 

FIGURE  1 


Unfortunately  the  specific  set  of  rules  for  E  to  maneuver 
by  to  achieve  the  survivability  depicted  in  Fig.  1  is  not 
known . 


II.   BACKGROUND 

A.   DISCRETE  EVASIVE  GAMES 

It  was  felt  that  a  two  dimensional  continuous  game  was 
too  difficult  to  solve.   So  the  approach  was  to  solve  a 
similar  but  easier  game.   It  was  assumed  that  the  ocean  was 
a  linear  set  of  discrete  points  and  that  E's  mobility  con- 
sisted of  being  able  to  jump  either  right  or  left  to  an 
adjacent  point  [Ref s ,  2  and  3].   Time-late  became  an  integer 
number  of  jumps  that  E  could  make  between  the  attack  de- 
cision and  warhead  detonation.   As  in  the  exact  problem  E 
had  no  knowledge  of  the  decision  to  attack. 

Discrete  problems  of  this  type  are  called  evasive  games 
and  are  classified  by  the  number  of  jumps  constituting  time- 
late.   A  one  step  discrete  evasive  game  means  E  can  move 
either  left  or  right  one  jump  prior  to  detonation.   Game 
theory  provides  an  immediate  solution  to  that  game;  the 
value  is  pj,=0,5.   Also  there  is  an  optimal  strategy  for  E 
in  that  game.   At  each  jump  E  should  go  left  with  probability 
one-half  or  right  with  probability  one-half.   Employment  of 
such  a  rule  for  each  jump  guarantees  E  a  survival  probability 
of  at  least  the  game  value,  regardless  of  P's  firing  rule. 

Having  easily  solved  the  one  step  game  interest  was 
focused  on  the  two  step  game.   In  this  game  E  was  allowed 
two  jumps  during  time-late.   The  solution  has  been  obtained, 
but  it  was  not  as  easy  to  achieve  as  that  of  the  one  step 
game  [Ref .2] . 


The  direction  of  the  analysis  was  clear.   Knowledge 
gained  from  solving  the  more  simple  games  would  be  a  step- 
ping stone  to  the  solution  of  the  more  difficult  games. 
Eventually  the  assumptions  of  linearity  and  discreteness 
could  be  relaxed.   Progress,  however,  in  solving  the  dis- 
crete games  has  been  very  slow.   Researchers  are  presently 
embroiled  in  solving  the  three  step  game. 

B.   A  CONTINUOUS  EVASIVE  GAME 

Washburn  in  Ref.  1  presented  a  different  approach  to 
analyzing  the  ICBM  versus  aircraft  carrier  problem.   He  de- 
veloped the  probability  density  function  of  a  particle 
moving  continuously  in  a  two  dimensional  medium  subject  to 
a  specific  set  of  maneuvering  rules .   That  approach  was 
distinct  from  previous  work  because  it  addressed  the  exact 
problem.   The  results  presented  in  this  paper  are  an  ex- 
tension of  that  approach. 

Washburn  proposed  the  following  strategy  for  E  and 
analyzed  the  evader's  subsequent  survivability,  first  as  a 
function  of  the  attacker's  strength  and  secondly  as  a 
function  of  a  specific  strategy  parameter.   E  was  to  select 
courses  from  a  uniform  distribution.   Course  changes  were 
to  occur  as  a  Poisson  process.   This  meant  the  time  between 
course  changes  would  be  exponentially  distributed. 

The  exponential  time  between  course  changes  has  an 
intuitive  appeal  because  of  the  memoryless  property  of  that 
distribution.   The  memoryless  property  is:  the  probability 
that  E  will  not  change  course  in  the  future,  given  that  he 
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has  not  changed  course  for  some  observation  period,  is  in- 
dependent of  the  length  of  that  period.   Therefore,  use  of 
the  exponential  time  distribution  should  serve  to  confound  a 
prediction  of  future  position  by  extrapolation,  regardless 
of  the  course  distribution  used.   For  example,  suppose  P's 
attack  decision  rule  required  that  he  observe  E  maintain  a 
constant  course  for  at  least  five  hours  prior  to  initiating 
an  attack.   Such  a  procedure  would  not  improve  P's  prob- 
ability of  a  kill  because  the  only  information  of  any  benefit 
at  the  time  an  attack  is  initiated  is  E's  position  and  last 
course. 

Figure  2  is  a  graph  of  the  evader's  maximum  attainable 
survivability  when  using  the  strategy  of  uniform  courses 

and  exponential  times  versus  P's  strength  for  the  optimal 

3 
exponential  parameter,  X.    For  comparison  the  upper  bound 

on  survivability  is  also  presented. 

"If  we  assume  the  parameters  A,v,TL  are  known  to  both 
sides,  then  the  evader  can  select  X   to  maximize  the  sur- 
vivability.  The  evader  will  clearly  be  in  trouble  if  he 
makes  X  too  small,  because  the  kill  probability  is  at  least 
exp(-X*TL).    On  the  other  hand  making  course  changes  too 
frequently  will  lead  to  a  density  function  that  is  highly 
peaked  at  the  origin,  which  is  equally  undersirable. . . " 
[Ref.  1], 
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The  exponential  parameter,  X,  is  the  inverse  of  the  mean 

time  between  course  changes. 

4- 
Exp(-X'TL)  is  the  probability  that  E  will  make  no  course 

change  during  TL. 


UPPER  BOUND 


0.2    0.4 


0.6    0.8 
A/tt  (vTL) 
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Minimum  pj^  for  Optimal  Parameter  A  Versus  Attacker 
Strength;  Strategy  Is  Uniform  Courses 
and  exp(A)  Times 

FIGURE  2 

"It  is  not  known  whether  or  not  the  optimal  strategy  for 
E  is  a  Poisson  strategy  of  the  type  just  considered,  or  even 
whether  the  uniform  distribution  on  angles  is  optimal  within 
the  class."   [Ref.  2] 
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III.   RESULTS  OF  SIMULATING  THE  CONTINUOUS  EVASIVE  GATIE 

The  goal  of  this  research  was  to  determine  if  the  uni- 
form course  rule  was  optimal  within  the  class  of  Poisson 
strategies  or  to  find  a  better  rule  if  it  was  not.   For 
another  rule  to  be  better  than  the  uniform  one  it  would 
have  to  enable  E  to  attain  a  higher  survivability  than  that 
shown  in  Fig.  2. 

The  methodology  of  solving  discrete  games  had  not  ad- 
vanced sufficiently  to  achieve  those  goals .   The  methods 
used  by  Washburn  were  also  not  suitable  because  of  the 
mathematical  difficulty  of  the  problem. 

The  Poisson  class  of  strategies  is,  however,  uniquely 
suited  for  analysis  by  computer  simulation.   This  is  due  to 
the  memory less  property  of  the  exponential  times  between 
course  changes.   The  results  of  such  a  simulation  are  pre- 
sented in  this  paper  in  the  form  of  improved  survivability 
for  the  evader. 

Five  course  change  rules  were  evaluated  in  the  simula- 
tion.  The  probability  density  function  of  each  is  pre- 
sented on  subsequent  pages  along  with  a  graph  of  the  re- 
sulting attacker's  pj^  versus  attacker  strength,  denoted  as 

o  5 
S=A/7T  (vTL) '^ .    These  graphs  are  for  representative  values 

of  the  parameter  a.   Alpha,  a,  is  the  product  of  time-late 

and  the  exponential  parameter  A. 


5 
Course  changes  were  assumed  to  be  independent  of  the 

underlying  Poisson  process. 
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The  simulation  established  a  nev7  lov/er  bound  on  the  pj, 
the  evader  could  yield  to  the  attacker.   Figure  3  is  a  graph 
illustrating  the  improvement  of  the  nev;  bound  over  that 
found  by  Washburn,   This  new  bound  is  the  least  lower  bound 
of  all  the  pj^  versus  S  curves  for  the  four  values  of  a   and 
all  rules  simulated.   Examples  of  such  curves  are  shown  in 
Figs,  5-9,   The  old  bound  in  Fig,  3  is  the  lower  bound  of 
the  pj,  curves  simulated  using  the  uniform  rule  exclusively. 

The  minimal  p,^  as  a  function  of  a,  over  the  five  rules, 
is  shov7n  in  Fig,  4  for  three  specific  attacker  strength 
levels.   The  curves  in  Fig,  4  confirm  the  statement  in  Ref. 
2  that  the  evader  should  turn  most  often,  that  is  a  should 
be  highest,  against  the  weakest  opponent. 

Comparisons  of  the  pj,  curves  of  the  different  rules 
have  shown  that  the  uniform  rule  is  not  optimal  for  all 
attacker  strengths  and  a  values.   For  example,  the  reverse 
course  rule  (see  Fig.  7)  was  shown  to  be  better  than  the 
uniform  rule  for  a=l.   The  difference  between  these  two 
rules  is  illustrated  in  Fig.  10.   Against  a  weak  opponent 
(strength  ^,6)  the  evader  would  do  better  to  use  the  re- 
verse course  rule  rather  than  the  uniform  rule  for  a=l. 
This  improved  survivability  occurs  because  the  reverse 
course  rule  created  a  "flatter"  evader  position  density. 

None  of  the  five  course  rules  evaluated  produced  a  con- 
sistently smaller  pj,  for  all  values  of  attacker  strengths 
and  a.   However,  the  left-right  rule  v/as  dominated  by  the 
oiiher  four  in  all  cases.   Each  of  the  five  course  rules  was 
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simulated  at  six  different  levels  of  a  by  varying  ET .   The 
graphs  of  pi^  versus  A/it(vTL)^  for  each  simulation  are  pre- 
sented in  the  Computer  Output  section  along  with  tabulated 
Pj^'s  for  specific  attacker  strength  levels  as  a  function 
of  the  course  rule  and  the  a  used. 


0.4- 


New  Bound  on  pj^  That  the  Evader  Can  Yield 
as  a  Function  of  Attacker  Strength 


FIGURE  3 
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Density  Function  of  Uniform  Rule 
FIGURE  5A 


a=1.0 


FIGURE  5B 
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Density  Function  of  Modified  Left-Right  Rule 


FIGURE  6 A 
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FIGURE  6B 
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Density  Function  of  Reverse  Course  Rule 


FIGURE  7 A 
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FIGURE  7B 
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FIGURE    8A 
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FIGURE  9A 
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FIGURE  9B 
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A  Comparison  of  the  Tradeoffs  Between  the  Uniform 
and  Reverse  Course  Rules  As  a  Function  of  Attacker  Strength 


FIGURE  10 
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IV.   SIMULATION  MODEL 

Let  E's  position  at  time  t  be  at  X=0 ,  Y=0  of  a  cartesian 
coordinate  system.   Then  rotate  the  axis  until  the  positive 
abscissa  is  aligned  on  and  pointing  in  E ' s  last  direction 
of  travel.   Then  regardless  of  the  course  rule  used  by  E 
there  is  at  least  exp(-a)  probability  that  he  will  be  at 
X=vTL,  Y=0  at  time  t+TL.   This  is  the  probability  that  E 
does  not  change  course  during  TL.   So  the  attacker's  prob- 
ability of  killing  E  by  centering  a  small  portion  of  his 
lethal  area  at  (vTL,0)  is  exp(-a).   Then  add  to  exp(-a) 
the  integral  of  E's  density  function  over  that  part  of  the 
uncertainty  circle  where  the  rest  of  the  lethal  area  is 
targeted  to  determine  the  total  pj^.   The  total  amount  of 
lethal  area  allocated  in  this  manner  is  the  numerator,  A, 
of  the  attacker's  strength  function  S. 

A  computer  program  to  simulate  the  above  computation 
was  written  in  FORTRAN  and  run  on  the  IBM  36  0.   The  program 
consisted  of  two  phases,  first  the  play  of  a  strategy  and 
secondly  the  scoring  of  that  play. 

A.   SIMULATION  OF  A  STRATEGY 

To  create  a  single  play  the  evader's  track  was  simulated 
from  an  initial  position  at  time  t  to  the  resultant  position 
at  t+TL.  The  track  was  the  result  of  a  specific  maneuvering 
strategy  being  simulated.  A  strategy  was  made  up  of  two 
decision  rules.  The  first  rule  determined  the  times  be- 
tween course  changes  which  were  exponential  random  variables 
with  mean  ET=1/A.   This  rule  was  common  throughout  all  of 
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the  simulations,  although  the  parameter  ET  was  a  variable. 

The  second  rule  determined  the  magnitude  of  subsequent 

course  changes.   The  only  other  kinematic  restriction  was 

a  constant  speed,  v,  for  the  evader. 

To  start  a  play  datum  was  initialized  by  setting  t=0 , 

X(0)=0,  Y(0)=0  and  an  initial  course,  c^ ,  was  selected  from 

the  uniform  distribution.   Then  a  sequence  of  exponential 

times  [tj^]-_T  vT/mT\   were  generated  and  a  sequence  of  course 

changes,  [c£]._^  ^  ,^^ .  ,   were  generated  using  the  course 

change  rule.   The  evader's  position  at  t=TL  was 

N-1 
X(TL)  =  X   v*t^»  (Cos  (c^)-Cos(c^)  )  +  vTL'Cos  (c^)    1.1 
i=l 


N-1 
Y(TL)  =  X  v*t.  •  (Sin(c.  )-Sin(c^)  )  +  vTL«Sin  (Cj,)    1.2 
i=l     ^ 


B.   SCORING  THE  PLAY 

Once  E's  position  at  t=TL  was  determined  that  observa- 
tion was  scored.   To  score  a  play  the  coordinates  of  E's 
position,  X(TL)  ,  Y  (TL)  ,  had  to  be  transformed.   The  purpose 
of  the  transformation  was  to  make  the  observed  position  in- 
dependent of  the  particular  initial  cource  c^ .   The  trans- 
formation was  a  rotation  of  the  coordinate  axis  about  the 
datum  so  the  positive  abscissa  would  be  aligned  in  the 
direction  of  the  initial  course.   The  transformed  position 
was : 


The  number  of  elements  of  this  sequence,  N(TL) ,  is  a 
Poisson  random  variable. 
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X'  (TL)  =  X(TL)Cos(c  )  +  Y(TL)Sin(c  ) 
Y'  (TL)  =  Y(TL)Cos(c  )  -  X(TL)Sin(c  ) 


2.1 
2.2 


An  example  of  this  transformation  is  shown  in  Fig.  11. 

If  E  had  not  made  a  course  change  during  TL  then  his 
transformed  position  would  have  been 

X'  (TL)  =  vTL  3.1 

Y' (TL)  =  0  3.2 

7 
A  grid  system  of  square  cells  was  placed  over  the 

playing  area  of  Fig.  IIB  and  a  determination  was  made  as  to 


R?VDIUS    =    vTL 


ORIGIN:    Datum 

B:    X(TL) ,    Y(TL) 

C  :  Initial  Course 

SAMPLE  PLAY  PRIOR 
TO  TRANSFORMATION 


ORIGIN:  Datum 

B:  X' (TL) ,  Y' (TL) 

A:  (vTL,0) 

SAMPLE  PLAY  AFTER 
THE  TRANS FOmiAT ION 


FIGURE  llA 


FIGURE  IIB 


which  cell  X' (TL) ,  Y' (TL)  was  in.   Each  cell  of  the  grid 
had  an  associated  value  which  v/ould  represent  the  number  of 
times  a  play  resulted  in  an  observation  in  that  cell.   V?hen 


The  cumulative  cell  area  was  noirmalized  by  the  factor 
7T  (vTL)  2. 
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the  appropriate  cell  was  deterroined,  for  the  play  being 
scored,  that  value  was  incremented  by  one. 

A  simulation  run  was  composed  of  16  810  plays  and  scoring 
iterations  for  a  strategy  utilizing  a  specific  course  change 
rule.   The  input  variables  for  a  run,  besides  the  course 
rule  were  v,  TL  and  ET.   At  the  end  of  a  simulation  run  the 
grid  system  was  a  two  dimensional  histogram  of  E's  position. 
The  frequencies  in  the  histogram  were  then  ordered,  accumu- 
lated and  normalized  to  achieve  cumulative  cell  probabilities 
The  ordering  corresponded  to  the  conservative  assimiption  that 
P  could  divide  his  lethal  area  and  target  only  those  cells 
with  the  higher  probabilities.   These  probabilities  were 
then  plotted  against  the  cumlative  cell  area  they  repre- 

o 

sented.    This  graph  was  labeled  p,  versus  attacker  strength 
and  was  the  primary  output  of  the  program. 

Thirty  simulation  runs  were  made  to  investigate  five 
different  course  change  rules  and  six  different  values  of 
ET.   For  all  runs  the  following  constant  values  were  main- 
tained, v=5  knots  and  TL=2  hours.   The  graphs  from  those 
runs  are  included  in  the  Computer  Output  section.   Each 
course  change  rule  was  simulated  many  different  times,  the 
only  difference  between  runs  being  the  pseudo  random  num- 
ber generator  seeds.   This  was  done  to  check  for  the  vari- 
ability of  the  pj^  graph  for  that  rule.   In  all  such  runs 


8 

The  amount  of  area  in  a  cell  was  denoted  as  cell  size  and 

egual  to  (2*TL/41)2,   tl  and  v  were  held  constant  there- 
fore cell  size  was  always  approximately  1/4  nm^. 


24 


the  resultant  graphs  were  so  similar  that  any  difference 
was  not  distinguishable.   For  simulation  run  required  two 
minutes  and  forty- five  seconds  using  125  K  on  the  IBM  360, 
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V.   CONCLUSION 

It  was  determined  that  the  uniform  course  change  rule 
is  not  optimal  within  the  Poisson  class  of  strategies .   Com- 
parison of  various  rules,  see  Fig.  10,  showed  that  under 
certain  conditions  a  single  rule  such  as  the  reverse  course 
change  rule  is  better.   Also  the  new  lower  bound  on  pj, , 
taken  over  all  rules  and  values  of  a  simulated,  is  an  im- 
provement on  the  lower  bound  achieved  from  only  the  uniform 
rule. 

Certainly  not  all  the  possible  course  change  rules  were 
simulated.   The  rules  evaluated,  however,  were  representa- 
tive of  the  broad  class  of  possible  rules.   When  compared 
to  the  uniform  rule,  all  but  one  of  the  other  rules  showed 
that  E  could  improve  his  situation  if  he  knew  P's  strength 
by  selecting  the  better  rule  for  that  encounter.   The  one 
rule  that  was  consistently  dominated  was  the  left-right 
rule. 

It  remains  unknown  whether  or  not  the  optimal  strategy 
is  a  Poisson  strategy. 
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APPENDIX   A 
FLOVJ   DIAGRAM 


START 


IX  —    55555 


ipseudo-random  number    seed,    an   input   variable 

> 


I  V  --  5-0   I 


c 


TL  ~  2.0 


1  evader  speed  in  knots  \ 

itime-late  in  hours    | 


ET  —  1.0 


,mean  time  between  course  changes,  an  inpu 

{variable __^ 


D 


PP(41,41) 
PP(i,j)  —  0 
for  all  i,j 


, playing  area  and  histogram;   each  cell  is 

'a  square  with  dimension  (2*V*TL)/41 


INDEX(41,41) 
INDEX(i,j)  —  0 
for  all  i,j 


—  I  INDEX  is  a  book-keeping  record  of  PP  cell 
.feasibility 


3 


rJINDEX(i,j)  ~  0  j5^ 


Is  PP(i,j)  an  exterior  cell 
J 


No 


INDEX(i,j)  —  ije-^^ls  PP(i,j)  an  interior  cell    ] 

No 


> 

X  —   0.0 

y  —   0.0 
Clock  ~   0.0 
c  —   0.0 

\ 

/ 

Itr         tt  — n  o^i- 

_ 

© 
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G> 


c  —  c  + c 

x 


Generate  a  course 

change   c  by  the 

rule  being  sim- 

lated     

TR 


|tR  —  TL-Clock 


No 


5  Generate  cw»u(0,2   ) 


Theta  —  c 


Generate  t*^Exp(l/ET) 


T  ~  min(TR,t) 


I 


Clock  —  Clock  +  t 


X  —  X  +  (V*TL)*Cos(c) 
y  —  y  +  (V*TL)*Sin(c) 


Is  Clock  =  TL 


.Yes 


X  ~  x*Cos(Theta)+y-Sin (Theta) 
y  ~  y-Cos(Theta)-x*Sin(Theta) 


Determine  which  column  of  PP  the  x  co- 
ordinate  is  in,  e.g.  column  1 


Determine  which  row  of  PP  the  y  coordi- 
nate  is  in,  e.g.  row  m 


PP(m,l)  ~  PP(m,l)  +  1 


STOP   ♦- 


^ 


Sort  and  order  PP 


Sum  PP 


A  Grap 
"\_i^ 


ph:   p,  versus  attacker  strength,  \ 
P|^  versus  normalized  lethal  area  \ 
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Subroutine  Trktm(ET,r , t) 

V 

|t ET*ln(r) 

1     ^ 

1 Return 

i 

Main 

1  Subroutine  Corse(r,c) 

> 

f 

|pi  —  3.1415 

' 

F 

f  —  2*Pi 

' 

' 

|a  —  0 

' 

' 

b  ~  f 

< 

r 

delta  c  —  r*(b-a)  +  c 

\ 

' 

c  —  c  +  delta  c 

^__  \ 

' 

c  —  c,Mod(f) 

y 

r 

Return 

N 

r 

Main 

Subroutine  Tcorse(r,c) 


s  —  0.2 


^ 
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(E> 


No 


*  Pi  ~  3.1415 


f  ~  2*Pi 


Is  r  between  0  and  0.25  inclusive 


Yes 


a  --  s/2 


b  --  l/(2*Pi)  -  s'^-Pi/A 


d r 


e  —  (b*"2)-A''<a-^d 


delta  c  —  (-b+Sqrt(e))/(2"-a) 


<D 


0- 


No 


Is  r  greater  than  0.25  and  less  than  or 
equal  to  0.5 


Yes 


a s/2 


b  —  l/(2--^Pi)  +  3*s''^Pi/4 


d  ~  (-s*(PiA*2)/4)  -  r 


e  ~  (b**2)  -  A'-^a^d 


delta  c  ~  (-b  +  Sqrt  (e)  )  /  (2'^-a) 


<D 
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^ 


No 


i 


Is  r  greater  than  0.5  and  less  than  or 
equal  to  0.75 


Yes 


rprime  —  r  -  0.5 


a  —  s/2 


b  —  l/(2-^Pi)  -  s*Pi/4 


d 


rprime 


e  ~  (b**2)  -  4=^a*d 


delta  c  —  (-b  =  Sqrt(e))/(2*a) 


delta  c  —  delta  c  +  Pi 


^ 


^  rprime  —  r-0. 5 


a s/2 


b  ~  l/(2*Pi)  4-  S-^s'^-Pi/A 


d  —  (-s''^(PiA*2)/4)  -  rprime 


e  --  (b**2)  -  4*a*d 


delta  c  —  (-b  +  Sqrt (e) ) / (2*a) 


delta  c  —  delta  c  +  Pi 


XD 
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1 


c 

—   c  +  delta   c        | 

\f 

c 

—   c,Mod(f) 

J, 

Return        J 


Main 


Subroutine  Pcorse(r,c) 


Pi  ~  3.1415 


f  ~  2*Pi 


s  ~  l/(?±->-<*2) 


Is  r  greater  than  0.5 


No 


>  a  —  s/2 


Yes 


3  —  -s/2 


delta  c  —  Sqrt(a) 


b  ~  l/(2*Pi)  +  3*Pi'^-s/2 


d s*(Pi>^-*2)  -  r 


e  —  (b**2)  -  A'^a^d 


delta  c-^   —    (-b+Sqrt(e))/(2*a) 


delta  C2  ~  (-b-Sqrt (e) ) / (2*a) 


5 


t 


32 


® 


i 


Is  Pi  less  than  delta  Ci  and  is 
delta  c-j^  less  than  or  equal  to  f 


Yes 


No 


delta  c  —  delta  c- 


I>1<^ 


delta  c  —  delta  c. 


|c  —  c  +  delta  c 

i 

c  —  c,Mod(f) 

> 

' 

Return 

1 

f 

Main 

T 


Subroutine  Trncrs(r 

,c) 

^ 

f 

Pi  ~  3.1415 

> 

f 

f  ~  2*Pi 

^ 

' 

a  ~  Pi/6 

I 


b  —  ll*Pi/6 

.  ,  y 

f 

delta  c  —  r*(b-a)  +  a 

^ 

' 

c  —  c  +  delta  c 

A 

c  ~  c,Mod(f) 

J 


1  Return 


Main 
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No 

Subroutine  Lrcors(r,c) 

i 

Pi  —  3.1415 

i 

f  --  2*Pi 

i 

c  —  c 

+  Pi/2 

^ 

Is  r  greater  than  0.5 

\ — 

Yes 

c  —  c  +  3*Pi/2 

a 

c   —   c,Mod(f) 


Return 


Main 
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column 

minimum 


Rule 


column 
minimum 


COMPUTER   OUTPUT 
TL/ET 


Rule 

1.0 

2.0 

3.0 

A.O 

Uniform 

.729 

.580 

.501 

.486 

Modified 
Left-Right 

.737 

.580 

.501 

.486 

Reverse 
Course 

.690 

.572 

.548 

.564 

Left-Right 

.917 

.768 

.650 

.588 

Truncated 

Uniform 

.697 

.556 

.509 

.517 

.690 


.556 


.501 


p,  for  attacker  strength  equal  0.2 
Table  I 
TL/ET 
1.0         2.0         3.0 


.838 


.768 


.752 


Pi^  for  attacker  strength  equal  0.4 
Table  II 


.486 


4.0 


Uniform 

.854 

.784 

.752 

.760 

Modified 
Left-Right 

.878 

.784 

.752 

.768 

Reverse 
Course 

.838 

.791 

.799 

.831 

Left-Right 

.964 

.909 

.870 

.854 

Truncated 
Uniform 

.846 

.768 

.768 

.815 

.760 


35 


TL/ET 
Rule  1.0        2.0        3.0        4.0 


Uniform 

.933 

.909 

.901 

.917 

Modified 
Left-Right 

.948 

.909 

.901 

.917 

Reverse 
Course 

.933 

.917 

.925 

.948 

Left-Right 

.987 

.987 

.980 

.987 

Truncated 
Uniform 

.933 

.909 

.909 

.933 

column  .933        .909        .901        .917 

minimum 

Pj^  for  attacker  strength  equal  0.6 

Table  III 
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f 


Pk- 

1.0 


0.8  - 


0.6 


0.4  - 


0.2  - 


0.2 


0.4 


0.6 


0.8 


A/Tr(v.TL)' 


1.0 


p^  versus  Attacker  Strength 
Course  Change  Rule:  Uniform 


ET  =  2.0 
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tr 


Pk 

1.0 


0.8 


0.6 


0.4 


0.2 


0.2         0.4 


0.6 


pu  versus  Attacker  Strength 
Course  Change  Rule:  Uniform 


0.8         1.0 
A/7r(v.TL)2 

ET  =  1.33 
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Pk 

1.0 


0.8  - 


0.6 


0.4 


0.2  - 


A/TT  (vTL)' 


pj^  versus  Attacker  Strength 
Course  Change  Rule:  Uniform 


ET  =  1.0 
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Pk 

1.0 


0.8 


0.6 


0.4  • 


0.2 


0.2 


0.4 


0.6 


0.8        1.0 
A/7T(v-TL)2 


Pl^   versus  Attacker  Strength 
Course  Change  Rule:  Uniform 


ET  =  0.1 
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Pk 

1.0 


0.8  - 


0.6  - 


0.4  - 


0.2 


0.2 


0.4 


0.6 


Pl^  versus  Attacker  Strength 
Course  Change  Rule:  Uniform 


0.8         l.Q 
A/Tr(v.TL)2 

ET  =  0.66 


41 


Pk 

1.0 


0.8  - 


0.6  - 


0.4 


0.2  - 


A/T^(v.TL)' 


Pl^  versus  Attacker  Strength 
Course  Change  Rule:  Uniform 


ET  =  0.5 


42 


Pk 
1.0 


0.8 


0.6 


0.4  - 


0.2  - 


0.2 


0.4 


0.6 


Pl^  versus  Attacker  Strength 
Course  Change  Rule:  Modified  Left-Right 


0.8        1.0 
A/-K  (v.TL)2 

ET  =  2.0 


43 


Pk 

1.0 


0.8  - 


0.6 


0.4 


0.2  - 


0.2 


0.4 


0.6 


P|^  versus  Attacker  Strength 
Course  Change  Rule:  Modified  Left-Right 


0.8        1.0 
A/T\  (v.TL)2 


ET  =  1.33 


44 


Pk 

1.0 


0.8 


0.6  - 


0.4  - 


0.2  - 


0.2 


0.4 


0.6 


0.8 


A/i^(v.TL) 


1.0 


pj^  versus  Attacker  Strength 
Course  Change  Rule:  Modified  Left-Right 


ET  =  1.0 


45 


Pk 

1.0 


0.8 


0.6  - 


0.4  - 


0.2  -1 


0.2 


0.4 


0.6 


0.8 


A/tC  (v.TL)' 


1.0 


Pj^  versus  Attacker  Strength 
Course  Change  Rule:  Modified  Left-Right 


ET  =  0.8 
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Pk 

1.0 


0.8  - 


0.6  - 


0.4  - 


0.2  - 


0.2 


0.4 


0.6 


Pj^  versus  Attacker  Strength 
Course  Change  Rule:   Modified  Left-Right 


0.8 


A/tt  (v-TL)' 


ET  =  0.66 


1.0 
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Pk 

1.0 


0.8  - 


0.6  - 


0.4 


0.2  - 


0.2        0.4 


0.6 


0.8 


A/Tr(vTL)' 


1.0 


p,  versus  Attacker  Strength 
^k 

Course  Change  Rule:  Modified  Left-Right 


ET  =  0.5 


48 


0.2 


0.4 


0.6 


0.8        1.0 
A/7^(vTL)^ 


Tpy^   versus  Attacker  Strength 
Course  Change  Rule:  Reverse  Course 


ET  =  2.0 


49 


0.2 


0.4        0.6 


0.8 


A/-K  (vTL) 


1.0 


p,  versus  Attacker  Strength 
Course  Change  Rule:  Reverse  Course 


ET  =  1.33 
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Pk 

1.0 


0.8  - 


0.6  - 


0.4  - 


0.2  - 


0.2 


0.4 


0.6 


P]f^  versus  Attacker  Strength 
Course  Change  Rule:  Reverse  Course 


0.8 


A/Tv  (v-TL) 


ET  =  1.0 


1.0 


51 


Pk 

1.0 


0.8 


0.6   - 


0.4 


0.2   - 


0.2 


0.4 


0.6 


0.8 


k/-K   (vTL) 


1.0 


pi^  versus  Attacker  Strength 
Course  Change  Rule:  Reverse  Course 


ET  =  0.8 


52 


Pk 

1.0 


0.8  - 


0.6  - 


0.4  - 


0.2  - 


0.2 


0.4 


0.6 


0.8 


k/K  (vTL) 


1.0 


Pi^  versus  Attacker  Strength 
Course  Change  Rule:  Reverse  Course 


ET  =  0.66 


53 


Pk 

1.0 


0.8  - 


0.6  - 


0.4  - 


0.2 


0.2        0.4        0.6        0.8 


p^  versus  Attacker  Strength 
Course  Change  Rule:  Reverse  Course 


1.0 


A/l^  (v-TL) 


ET  =  0.5 


54 


Pk 

1.0 


0.8 


0.6 


0.4  - 


0.2 


0.2 


0.4 


0.6 


Pj^  versus  Attacker  Strength 
Course  Change  Rule:  Left-Right 


0.8 


A/ 7%  (vTL) 


ET  =  2.0 


1.0 
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Pk 

1.0 


0.8 


0.6 


0.4  ■ 


0.2  ■ 


0.2        0.4 


0.6 


0.8         1.0 
A/ti(vTL) 


p^  versus  Attacker  Strength 
Course  Change  Rule:  Left-Right 


ET  =  1.33 


56 


Pk 

1.0 


0.8  - 


0.6  - 


0.4 


0.2  - 


0.2 


0.4 


0.6 


0.8 


A/ 71  (vTL)' 


1.0 


Pj^  versus  Attacker  Strength 
Course  Change  Rule:  Left-Right 


ET  =  1.0 


57 


Pk 

1.0 


0.8  - 


0.6  - 


0.4  - 


0.2  - 


0.2        0.4 


0.6 


P|^  versus  Attacker  Strength 
Course  Change  Rule:  Left-Right 


0.8        1.0 

k/li   (v.TL)2 

ET  =  0.8 
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Pk 

1.0 


0.8 


0.6 


0.4 


0.2 


0.2 


0.4 


0.6 


Pj^  versus  Attacker  Strength 
Course  Change  Rule:  Left-Right 


0.8 


A/-K  Cv«TL)' 


ET  =  0.66 


1.0 


59 


Pk 

1.0 


0.8 


0.6 


0.4  - 


0.2 


0.2 


0.4 


0.6 


0.8 


1.0 


A/Ti  (vTL) 


p,  versus  Attacker  Strength 
k 


Course  Change  Rule:  Left-Right 


ET  =  0.5 


60 


0.2 


0.4 


0.6 


0.8 


A/n  (v-TL) 


1.0 


Pj^  versus  Attacker  Strength 
Course  Change  Rule:  Truncated  Uniform 


ET  =  2.0 


61 


u. 


1.0 


0.8- 


0.6  - 


0.4  - 


0.2  - 


A/-K  (v-TL) 


p^  versus  Attacker  Strength 
Course  Change  Rule:  Truncated  Uniform 


ET  =  1.33 
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0.4  - 


0.2  ' 


0.2 


0.4 


0.6 


Pj^  versus  Attacker  Strength 
Course  Change  Rule:  Truncated  Uniform 


0.8         1.0 
A/7\  (vTL)2 

ET  =  1.0 


63 


0.4  - 


0.2  - 


0.2 


0.4 


0.6 


0.8 


A/-R  (vTL)' 


1.0 


P]^  versus  Attacker  Strength 
Course  Change  Rule:  Truncated  Uniform 


ET  =  0.8 


64 


0.2 


0.4 


0.6 


0.8 


A/tc  (vTL) 


1.0 


p,  versus  Attacker  Strength 
Course  Change  Rule:  Truncated  Uniform 


ET  =  0.66 


65 


Pk 

1.0 


0.8  - 


0.6  - 


0.4 


0.2  - 


0.2 


0.4 


0.6 


0.8         1.0 
A/7T(v.TL)^ 


Pj^  versus  Attacker  Strength 
Course  Change  Rule:  Truncated  Uniform 


ET  =  0.5 
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COMPUTER   PROGRAM 

A   SIMULATION   OF   THE    CONTINUOUS    EVASIVE 
GAME    FOR  THE    POISSON    CLASS    OF    STRATEGIES 

INTEGER    *    2    PP, INDEX 

INTEGER   *    2    CUMSUM 

INTEGER    *    2    NCOUNT 

INTEGER    *2    XP,QQ 

INTEGER   *    4    GRIT'S 

INTEGER    *    A    TOTCEL 

REAL    *    A    IN 

REAL    LABEL/'  '/ 

REALMS    ITITL^(12)/'CL0THIER« .'       THOMAS'.'       J,  ',3*' 

X  ' ,'STRATEGY','       TRNCRS ' , '       (LAMDA','    *    TL )=     ', 

X2*»  4.0*/ 

C  BY    MAKING   THE    HISTOGRAM    INTEGER    STAR    TWO    THE    AMOUNT 

C  OF    STORAGE    REQUIRED    FOR    THE    PROGRAM    WILL     BE    REDUCED. 

C 
C 

DIMENSTGN    PP(41,41 ) , INDEX ( 41 , 41 ) 

DIMENSION    XP(1681) ,00(41 ,41) 

EQUIVALENCE  (PP (  1 , 1  ) , XP( 1 ) ) 

DIMENSION  EEX(2) , EEY (2 ) , EX (692 ) , EY( 692 ) » EXX( 692 ) 

DIMENSION  EYY(  692) 
C 
C 

•C  PP    IS    THE    PLAYING    AREA    AND    HISTOGRAM,       EACH    CELL    IN 

C  PP    REPRESENTS    AN    AREA    OF    SIZE    CELSIZ.       INDEX    IS    A 

C  BOOK-KEEPING    ARRAY    USED    TO    DETERMINE    WHICH    CELLS    IN 

C  PP    ARE    FEASIBLE,    E.G.    IF    IND£X(I,J)=0    TH!HN    PP(I,J)    IS 

C  AN    INFEASIBLE    CELL.       SIMILARLY    IF    INDEX(I,J)=1    OR    2 

C  THEN    PP(I,J)     IS    A    FEASIBLE    CELL. 

C 

C  PP    IS    A    SQUARE    WITH    THE    TIME-LATE    CIRCLE    INSCRIBED    IN 

C  IT.       THEREFORE    PP    HAS    SOME    CELLS    THAT    ARE    NOT    INSID^: 

C  THE    TIME-LATE    CIRCLE    AND    FOR    THOSE    CELLS    TH^IR    INDEX 

C  VALUE    IS    ZERO.       FOR    THOSE    CELLS    IN    PP    THAT    ARE    CON- 

C  TAINEO    ENTIRELY- IN    THE    CIRCLE    THE    ASSOCIATED    INDEX 

C  VALUE    IS    ONE.       AND    LASTLY    THERE    ARE    THE    CELLS    IN    PP 

C  THAT    ARE    ON    THE    EDGE    OF    THE    CIRCLE,    THESE    ARE    CALLED 

C  EDGE    CriLLS   AND    HAVE    THE    INDEX    VALUE    OF    TW3. 

C 

c 

NCCUNT=0 

PI=3. 1415926 

DATA    IX,V,TL,ET/37915,05.0,02.00,00.5/ 
C  VIS    PARTIClI    SPEED 

C  TL     IS    THE    TIME-LATE 

C  ET    IS    THE    MEAN    TIME    BETWEEN    COURSE    CHANGES 

C  IX    IS    THE    SEED    FOR    THE    PSEUDO-RANDOM    NUMBER    GENERATORS 

C 

c 

DELTAX=(2.C^V*TL)/41.0 
C  DELTAX    IS    THE    HORIZONTAL    DIMENSION    OF    A    CELL. 

C 

CELSIZ=DELTAX*t2 
C  CELSIZ    IS    THE    AMOUNT    OF    REAL    AREA    CONTAINED    IN   A    CELL. 

C 

K  =  20 

C  THE    NUMBER    OF    CELLS    IN    EITHER    DIMENSION    IS    ALWAYS    A 

C  CONSTANT    41.       THIS    CONSTANT     IS    EXPRESSED    AS    AN    ODD 

C  INTEGER    AND    WRITTEN    IN    THE    FOLLOWING    MANNER, 

C      Al  =  2  *  (K)  +  1.   THEREFORE  K  =  20. 

C 

C 

C      ZERO  OUT  THE  HISTOGRAM 

DO  5   1=1,41 

DC  4   J=l,41 

INDEX(  I, J  )  =  0 
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PP(I,J)=0 
00(1, J)=0 

4  CONTINUE 

5  CCNTINU2 
C 

c 

C  NOW    TO    CHECK    EACH    OF   THE    CELLS     IN    THE    HISTOGRAM    AND 

C  DETERMINE    ITS   ASSOCIATED    INDEX    VALUE. 

C 

c 

C  QUADRANT    I 

DO    42   J=21,41 
DO   41    1=1,21 

0UT=(V=*^TL-J-^DELTAX)+*2+(  ( I -1 )  *  DELTAX-V*TL)  **2 
IN  =(V*TL-(J-1  )=^DELTAX)v*2+(  I  =^D£LTAX-V=pTL  )  **2 
RADS0=(V=^TL)=^-*^2 

IF( (OUT.LE.RADSO).ANO. (IN. LE.PADSO) ) I NDEX ( I , J )=1 
.IF(  (OUT.GT.PADSO).  AND.dN.GT.PADSQ  ))IND.:X(  I,J)=0 
IF( (OUT.GT.RADSO).ANO.(IN.LE.RADSO) )IND;X(  I, J)  =  2 

41  CONTINUE 

42  CONTINUE 
C 

c 
c 

C      QUADRANT  II  CHECK 
DO  44  J =1,21 
DO  43  1=1,21 

OUT=(V=^TL-(  I-l  )=^DELTAX)*-*2  +  (V*TL-(J-l)*DELTAX)**2 
IN  =  (  V^TL-r^D'':LTAX)**2+(V=^TL-J=^^DELTAX)*=p2 
RADSQ=(V*TL  )-^*2 

IF( (QUT.LE.RADSO). AND.(IN.LE.RADSO) ) INDEX ( I, J)=l 
IF(  (nUT.GT.RADS0).AND.(IN.GT.RADS0))IND:3X(  I  .  J  )  =0 
IF( (OUT.GT.RADSQ).AND.(IN.LE.RADSQ) ) INDEX ( I,J)=2 

43  CONTINU" 

44  CONTINUE 
C 

c 

C  QUADRANT  III 
DO  46  J=l,21 
DC  45  1=21,41 

OUT=(V=^TL-(  J-1)*DELTAX)**2+(I*DELTAX-V*TL)**2 
IN=(V*TL-J*0ELTAX)**2+( ( I -1 )*DELTAX-V*TL )**2 
RADSQ=(  V^TL)-?^2 

IF( (OUT. LE. RAD  SO ).AND.( IN.LE.RADSO)  )INO£X(I , J)=l 
IF( (OUT.GT.RADSQ).ANn.(IN.GT.RADSQ) )INDEX(!,J)=0 
IF( (OUT.GT.RADSO).AND. (I  N.  LE.RaDSO) )INDEX( I , J )  =2 

45  CONTINUE 

46  CONTINUE 
C 

c 

C      QUADRANT  IV  CHECK 
DO  48  J=21,41 
DO  47  1=21,41 

OUT=(V^TL-I"  DELTAX  )^*2  +  (  V*TL-J*OELT  AX  )*'^2 
IN  =(  V»TL-(  I-l)-^DELTAX)4^2+(  V'feTL-(  J-l)*DeLTAX)-i'i^2 
RADSQ=(V-^TL)*--2 

IF( (OUT. LE.RADSO). AND. ( I N. LS .RADSO ) ) INDEX ( I , J )=1 
IF( (OUT.GT.RADSQ).AND. (IN.GT.RADSO) )INDSX(I , J ) =0 
IF( (nuT.GT.RAOSO).AND. ( IN.LE.RADSO) ) INDEX (  I, J  )  =  2 
A7  CONTINUE 
48  CONTINUE 

C 

c 

C  NOW    TO    SIMULATE    THE    PARTICLE    MOTION    UNDER    THE 

C  SPECIFIED    RULES.       THIS    IS    THE    MAIN    DO    LOOP    IN    THE 

C  PROGRAM.       FIRST    A    TIME    UNTIL    THE    NEXT    COURSE    CHANGE 

C  IS    GENERATED.       THEN    A    COURSE    TO    BE    STEERRED    IS    GEN- 

C  ERATED.       XX    AND    YY    ARE    THE    TWO    COMPONENTS    OF    POSITION 

C  ADDED    DUE    TO    A    COURSE    AND    TI  m*^       SEGMENT.        X    AND    Y    ARE 

C  THE    UPDATED    POSITION    ^RQM    DATUM.       AT    TIME-LATE    THE 

C  UPDATED    POSITION    IS    CONVERTED    TO    A    CELL     POSITION    AND 

C  THAT   CELL    VALUE     IS    INCREMENTED    BY    ONE. 
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c 
c 

DO  3  J  =  l, 16810 

X=0.0 

Y=0,0 

CLOCK=0.0 

C=0.0 

TR=TL-CLOCK 

CALL  RANDU(IX,IY,R) 

IX  =  IY 

CALL  CORSE(R,C) 

THETA=C 

GO  TO  66 
1  TR=TL-CLOCK 

CALL  RANDU(IX,IY,R) 

IX  =  IY 

CALL  TPNCRS(R,C) 
66  CALL  RANDU(IX,IY,R) 

IX=IY 

CALL  TRKTM(?T,R,T) 

TT=AMIN1(T,TR) 

CALL  MOV^(V,TT,C,XX,YY) 

CLOCK=CLOCK+T 

X=X+XX 

Y=Y+YY 

IF(CLOCK.LT.TL)';0    TC    1 

XPRIV|^=Xi=CnS(THETA)+Y^SIN(  TH?TA) 

YPRIM"=Y*COS{THnTA)-X*SIN(THHTA) 

X=XPRIME 

Y=YPRIM? 

CALL    RANOUdXt  IY,R) 

IX=IY 

C0RR=1.0 

CALL    HISTO(K,DELTAX,R,CORRtX,I ) 

L  =  I 

CALL    RANDU( IX, IY,R) 

IX  =  IY 

CGRR=-1.0 

CALL    HISTG(K,DELTAX,R,CORRTYtI ) 

N  =  I 
C  THH    PARTICULAR    CELL    I,  J    IS    DET^RMIN^D    BY    TWC       SUCCESS- 

C  IV£    CALLS    OF    HISTO.       IT    IS    POSSIBLE    A    POSITION    X,Y 

C  CCULD    BE    DSTERMIM^D   TO    Be    IN    AN    INFEASIBLE    CELL. 

C  IF    THIS    HAPPENS    TH5M    RESOLV    WILL    BE    CALLED    TO    CHANGE 

C  THAT    CELL    ASSIGNMENT    IN    A    PRcSCRIBEO    MANNER    TO    A 

C  FFASIBLS    CELL. 

IF(INDEX(M.L ).EO.O)GO    TO    31 

30  PP(M,L)=PP(M,L )+l 
GO    TO    3 

31  CALL    RESOLV(L,M, I,N) 
NC0UNT=NC0UNT+1 

L  =  I 
M=N 

GO   TO    30 
3    CONTINUE 
C 

NUPC5L=0 

NUPEDC=0 

CRITE2=0.0 

DO    52    1=1,^1 

DO    51    J=1t41 

IF(PP(  I»  J)  .EO.OGO    TO    51 

NUPCEL=NUPCEL+1 

IFdNDEXd,  J)  .NE.2)G0    TO    51 

NUPE0C  =  N')PEDC+1 

CRITE2=CRIT-2+PP(I t J) /I 6 81 0.0 

51  CONTINUE 

52  CONTINUE 
CRITE3=NUP^DC 
RNUPED=NUP:-DC 
RNUPCE  =  NUP'^^-L 
CRITE1=RNUP£D/RNUPC5 
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NUPINC=NUPCFL-NUPEOC 
212    WR!T'=(6,203) 

203  FCRMATC • 1' t60X,' PR       PAGE    1«) 
DO    205    1=1.41 
WRITE(6,204)     ( PP ( I , J ) » J=l t 15 ) 

204  FORMAT    (  •     '  ,1517) 

205  CONTINUE 
WRITE  (6,206) 

206  FORMAT(  'l' ,60X,'PP   PAGE  2M 
00  208  1=1,41 
WRITE(6,207) (PP(I, J),J=15,29) 

207  FORMAT( •  ',151 7) 

208  CONTINUE 
WRIT^     (6,209) 

209  FORMAT(  'l*  ,60X,  "PP       PAGE    3M 
DO    211    1=1,^1 
WRITE(6,210) (PP( I, J) ,J=29, 41) 

210  FORMAT( •     ',1317) 

211  CONTINUE 
C 

C      NOW  THAT  THE  HISTOGRAM  IS  BUILT  ITS  ELEMENTS   ARE 

C      ORDERED  SO  THE  HIGHEST  VALUES  ARE  IN  THE  LOWEST  NUM- 

C      BERED  CELLS. 

C 

r 

''lOCO  IC=0 

DO  1010  1=1,1680 

IFCXPd)  .GE.XP(I+1  )  )  GO  TO  1010 

MAX=XP( I+l) 

XP(I+1)=XP( I) 

XP(I)=MAX 

IC  =  10 
1010  CONTINUE 

IF(IC.NE.O)GO  TO  1000 

DO  2000  1=1.41 

DO  2010  J=l,41 

OQ(I  ,J) =PP( J,I  ) 
2010  CONTINUE 
2000  CONTINUE 

DC  1030  1=1,41 

DO  1020  J=l,41  . 

PP(I ,J)=00(I,J  ) 
1020    CONTINUE 
1030    CONTINUE 
C 

c 

C  AFTER    ORDERING    THE    CELL    VALUES    THEY    ARE    ACCUMLATED    SO 

C  EACH    CELL    CONTAINS    THE    SUM    OF    ITS    OWN    VALUE    AND    ALL 

C  THE    VALUES    PRECEDING    IT. 
C 

c 


c 
c 


CUMSUM=0 

DO    11    1=1,41 

DO    10    J=l,41 

CUMSUM=CUMSUM+PP(I  , J) 

PP( I , J)=CUMSUM 

10  CONTINUE 

11  CONTINUE 


NUMEXC=0 

NUMINC=0 

NUMEDC=0 

DO   33    1=1 ,41 

DO    32    J=l,41 

IF    (I NO EX (  I, J) .EQ.0)NUMEXC  =  NUMEXC  +  1 

IP    (IND^Xd  ,J)  .EQ.l  )NUMINC  =  NUMINC  +  1 

IF    (  INDEX(  I ,  J  )  .£0.  2)  NUMEDC  =  NUMEDC-H 

32  CONTINUE 

33  CONTINUE 
NUMC?L=NUM!NC+NUMEDC 
TOTCEL=NUMEXC+NUMEDC+NUMINC 
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c 
c 


c 
c 


c 
c 

C  NOW    TO    LOAD    THE  FIRST    900   CELL    VALUES    INTO    THE    PLOT- 

C  TING    VECTOR    EY.  THE    CUMULATIVE    NUMBER    OF    LOOKS    OR 

C  KILL    POWER    WILL  BE    PLOTTED    ON    THE    X    AXIS. 
C 

c 

EEX( 1)=0.0 
EEX(2)=1.0 
EEY{1)=0.0 
EEY(2)=1.0 

DO    15    1=1,16 

DO    14    J=l,41 

N  =  (l-l)^41  +  J 

EY(N)=PP(I»J)/16810.0 

EX(N)=N/1385.0 

14  CONTINUE 

15  CONTINUE 
DO    16    J=l,36 
1  =  17 

N=(I-1)*41+J 
EY(N)=PP(I ,J)/16810.0 
EX(N)=N/1385,0 

16  CONTINUE 

DO    17    J=37,41 

1=17 

N=( I-l)^41+J 

L=N-692 

EYY(L)=PP{ I, J)/16810.0 

EXX(L)=N/1385.0 

17  CONTINUE 
DO  19  1=18,33 
DO  18  J=1,A1 
N=(T-1)-^41+J 
L=N-692 

EYY(L)=PP(I,J)/16810.0 
EXX(L)=N/1385.0 

18  CONTINUE 

19  CONTINUE 
DC    20    J=l,31 
1=34 

N=(I-1)*41+J 
L=N-692 

EYY(L)=PP( I, J) /16810.0 
EXX(L)=N/1385.0 

20  CONTINUE 
WRITE(6,606)    NCOUNT 

606  F0RMAT(//,25X, 'NCOUNT  'OUALS   ',116) 

WRITE    (6,100)     NUM.1XC,NUMEDC,NUMINC,T0TCEL 

100  FORMAT  (/, lOX, •NUMEXC=', 116, 5X, 'NUMEDC=' ,1I6,5X,« NUMIN 
XC  =  '  ,1T6,5X,'T0TCEL=' ,116) 

WRITE (6, 101)    CRITE3,CRITE1.CRITE2 

101  FORMAT  (//,10X, • NUMBSR  OF  POSITIVE  EDGE  CELLS  IS  ',11 
X6,//,10X,'THc  PERCENTAGE  OF  POSITIV:^  CELLS  WHICH  ARE  E 
XDGE  CELLS  IS  • , IF  4. 3. // , lOX , • THE  AMOUNT  OF  PRQBABILIT 
XY    IN   THE    EDGE    CELLS    IS       ',1F4.3) 

CALL  DRAW  ( 2 ,EE X  ,EEY ,1 ,0 , LABEL , ITITLE ,0 ,0 ,0 ,0 ,0 , 0 , 9, 9, 
X0»LAST) 

CALL  DRAW  (692 , EX, EY ,2 ,0 , LABEL , IT ITLE, 0, 0, 0, 0, 0, 0, 9, 9, 
XCLAST) 

CALL  DRAW  (692,EXX,EYY,3,0,LABEL,ITITLE,0,0,0,0,0,0,9, 
X9,0,LAST) 

STOP 

END 
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StBROUTINE    RANDUd  XtlY.YFL) 

IY=IX*65539 

IF(IY) 5,6,6 

5  IY=IY+2147483647+1 

6  YFL=IY 
YFL=YFL*.4656613E-9 
RETURN 

END 


SUBROUTINE  TRKTM (E T ,R, T) 

RA=1.C-R 

T1=-ST*AL0G(R) 

T2=-ST:^AL0G(RA) 

T=T1 

RETURN 

END 


SUBROUTI N^    MOVH( VSL,TIME,ANG,X,Y) 

X=Vf^L^TTME*COS(ANG) 

Y=VEL*TIME*SIN(ANG) 

RETURN 

END 


SLBROUTINE    HISTO(K,DSLTAX,R,CORR,X,I ) 

UP3ND=(  {  2*K  +  l);^cr:LTAX)/2.0 

L  =  2*K-H 

XX=(X+(C0RR='^UP8ND)  )/DELTAX 

XX=CORR-^=XX 

Xl=AMaD(XX,l .0 ) 

M=INT( XX) 

IF(M.eO.O)GO    TG    1 

IF(M.EQ.L)GO    TT    2 

IF( X1.N5.0.0)G0    TO    3 

1F(R.GT.0.5)G0    TO    3 

I=M 

GO    TO    4 

1  1  =  1 

GO    TG    4 

2  I  =  L 

GO   TO   4 

3  I=M+1 

4  RETURN 
END 


SUBROUTINE    RESOLVd  ,J,L.M) 

IF(!.GT,21)G0    TO    2 

IF( J.GT.21)G0    TO    1 

L=I  +  1 

M  =  J  +  1 

GO    TO    4 

L=I  +  1 

^'=J-l 

GO    TO    4 

IF( J.GT.21)G0    TO    3 

L  =  I-1 

y  =  J+l 

GO    TO    4 

L  =  I-1 

H  =  J-l 

RETURN 

END 
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SUBROUTINE    CORSE(RtC) 

PI=3. 1415926 

A  =  0.0 

B=2.0*PT 

F=2.0*P! 

DELTAC=R*(B-A)+A 

C=C+DcLTAC 

C=AMnD(C,F) 

RETURN 

END 

SUBPOUTINc    TCORSS(R,C) 

C  R    IS    THE    UNIFORM    VARIATE    USED    TO    GENERATE    THE    NEXT 

C  COURSE    CHANGE, 
C 

C  C    IS    THE    PRESENT    COURSE    AND    WILL    BE    RETURNED    AFTER 

C  ADDING    A    NSW    HEADING    CHANGc. 
C 

C  DELTAC    IS    THE    AMOUNT    OF    THE    NEW    HEADING    CHANGE 
C 

S=0,20 
C 

C  S    IS    THP   JE^PZ^.    SLOPE    AND    IS    BETWEEN    0   AND 

C  2/(PI*>^2)    =    0.292.       NOTE    »    S    CANNOT    EQUAL    0. 


C 

c 

c 

c 

c 

c 
c 


c 
c 


C 
C 


c 

c 


PT=3. 1415927 

IF(0.0.LE.R.AND.R.LE.0.25)G0  TO  1 
IF(0.25.LT.R.AND.R.LE.0.5)G0  TO  2 
IF(0.5.LT.R.AND.R.LE.0.75)G0  TO  3 
IF(0.75.LT.R.AND.R. LE.1.0)G0    TO    4 


A=S/2.0 

B=(1.0/{2.0*PI  )  >-(  S-*PI/4.0) 

D=-1.0^=R 

^  =  8=^^=^  2-4.0* A*D 

DELTAC=(-1.0'5'B+S0RT{E)  )/(  2.0-^A) 

GO    TO    5 


RPRIME=R-0.5 

A=S/2.0 

B=(1.0/(2.n*PI)  )-{  S'^PI/4.0) 

D=-1.0^<RPRIMc 

E  =  B**2-4.0*A^-D 

DELTAC=(-l.C-^8+SQRT(E)  )/(2.0^A) 

DELTAC=DPLTAC+PI 

GO    TO    5 


A=(-1.0''*S)  /2.0 

B=(1.0/(  Z.O^^PI  )  )+(  3.0*S*PI/(  4.0)  ) 

D=-1.0-M  (S':PI^<-=2  )/4.0  )-R 

E=B*:^:2-4.0*A*D 

DELTAC=(-1.0'»-3+SQRT(E)  )/(2.0-*A) 

GO    TO    5 


RPRIME=R-0.5 

A  =  (-1.0:*  S)  /2.0 

B=(1.0/(2.0-*^FI  ))+(  3.0*S*PI/(  4.0)  ) 

D=-l  .0^  (  (S^P!^^^^2  )/4.0  )-RPRIME 

E  =  B:*v2-4.0^A'*=D 
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DELTAC=(-1.0*B+S0RT(E) )/(2.0*A) 

DELTAC=DELTAC+PI 

GO    TO    5 


c 

5  C=C+DELTAC 

c 

F=2.0*PI 

c 

C=AMOD(C,F) 

c 

RETURN 

END 

SUBROUTINE    PCORSE(R,C) 
PI =3. 141 5926 
F=2.0*PI 
S=1.0/(o!**2) 
IF(R.GT.0.5)GG    TO    1 
A=2.0^R/S 
DELTAC=SQRT(A) 
GO    TO    3 

1  A=-1.0^'S/2.0 
B=(1.0/{2.0^-PI  )  )+(  3. 0*P  1=^5/2.0) 
D  =  -1.0'^S>MP!*^2)-R 
E=(e-v*2)-4.0*A*D 

deltc1  =  (-b  +  sq'.t(=)  )/(2.c*a) 
deltc2=(-b-s0p.t(s)  )  /(2.0^a) 
ip(pi.it.d5ltci.aNo.deltci.le.f)go  to  2 
celtac=d^ltc2 

GO    TO    3 

2  DELTAC=D1LTC1 

3  C=C+D'LTAC 
C=AMOD(C,F) 
RETURN 

END 


SUBROUTINE    LRCORS(R,C) 
PI=3. 1415926 
F  =  2.0^'  PI 

1F(R,GT.0.5)G0    TO    1 
C=C+(PT/2.0) 
GO    TO    2 

1  C=C+((3.0-^PI  )/2.0) 

2  C=AMaD(C,F) 
RETURN 

END 


SUBROUTINE    TRNCRS(R,C) 

PI=3. 1415926 

F=2.0^PI 

A  =  P!/6.0 

B=(11.0^PI )/6.0 

DELTAC  =  Rt:(B-A)+A 

C=C+DELTAC 

C=AMOD(C«F) 

RETURN 

END 
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are  independent  exponential  random  variables.   An  improvement  over 
a  previously  established  least  upper  bound  on  evader  survivability 
is  presented.   The  bulk  of  this  paper  is  the  presentation  of  evader 
survivability  as  a  result  of  the  maneuvering  strategy  he  employs. 
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